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Abstract 

We recompute the functions describing the collinear factorization of one-loop amplitudes using 
the unitarity-based method. We present the results in a form suitable for use as an ingredient in 
two- loop calculations. We also present a function summarizing the behavior at one loop in both 
the soft and collinear limits. 
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1. Introduction 



The properties of non-Abelian gauge-theory amplitudes in soft and collinear limits play an important 
role both in formal proofs and in explicit calculations. The singularities that arise in these limits create 
technical obstacles to calculations of amplitudes and cross sections. Their universal structure, however, 
provides a means of resolving these problems. The universality of the limits is furthermore essential to the 
applicability of quantum chromodynamics to short-distance physics. It also provides a practical means for 
checking calculations of higher-point amplitudes. Explicit calculations of the functions governing these limits 
thus play an important role in the program of precision quantum chromodynamical calculations. 

In these singular limits, on-shell amplitudes factorize into sums of lower-point amplitudes multiplied by 
universal functions. In the collinear limit, these functions are called splitting amplitudes. At tree level, one 
may derive the splitting amplitudes from a string representation [1] or from the Berends-Giele recurrence 
relations [2]. The squares of these tree-level functions, summed over helicities, yield the Altarelli-Parisi 
kernels [3,4]. 

Beyond leading order, the splitting amplitudes do not yield the Altarelli-Parisi kernel directly, and the 
relation between the two quantities remains to be clarified. The explicit forms of the splitting amplitudes 
at one-loop through O(e ) (in dimensional regularization, with D = 4 — 2e) have previously been extracted 
by taking the collinear limit explicitly in various five-point amplitudes [5,6], or from an analysis of one-loop 
integrals [7]. In addition, Bern, Del Duca, and Schmidt [8] have recently given an expression for the gluon 
splitting amplitude to all orders in the dimensional regulator e. 

Many of the complicated one-loop computations of recent years have been performed using the unitarity- 
based method developed by Bern, Dixon, Dunbar, and one of the authors [9,10]. The method also leads 
to a concise proof [11] of collinear factorization to all loop orders in gauge theories, and as a by-product 
provides a compact formula for computing the splitting amplitudes explicitly. In this paper, we use this 
formula to compute the various one-loop splitting amplitudes, thus providing complete forms to all orders 
in e for both gluon and fermion external states. The splitting amplitudes presented here will be useful in 
computations of physical observables such as jet cross-sections, at next-to-next-to-leading order. In addition 
to the one-loop splitting amplitudes, for such purposes one would also need the tree-level functions describing 
factorization when three particles become collinear, computed by Campbell and Glover [12] and by Catani 
and Grazzini [13]. 

The paper is organized as follows. We review the structure of collinear factorization in the next section. 
We then present a detailed example calculation, the gluon splitting amplitude with selected external helicities, 
in section 3. We present the complete set of gluon and fermion splitting amplitudes in section 4, and discuss 
various checks and comparisons in section 5. In section 6, we use these splitting amplitudes to obtain a 
function which interpolates smoothly between the collinear and soft limits. 
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2. The Structure of Collinear Limits 



The properties of gauge theories are easiest to discuss in the context of a color decomposition [14]. At 
tree level, for all-gluon amplitudes such a decomposition takes the form, 

AT e ({k i ,X i ,a i })=g n - 2 £ T^TMD ■ ■ • T<M-> )< re >(l A \ . . . ,n A ")) , (2.1) 

<res n /z n 

where S n /Z n is the group of non-cyclic permutations on n symbols, j s denotes the j-th momentum and 
helicity, and g is the gauge coupling. We use the normalization Tr(T a T b ) — 5 ab for the generators of SU(N). 
One can write analogous formulas for amplitudes with quark-antiquark pairs or uncolored external lines. The 
color-ordered or partial amplitude A n is gauge invariant. 




Figure 1. A schematic depiction of the collinear factorization of tree- level amplitudes, with the 
amplitudes labelled clockwise. 

In the collinear limit, k a | fc& of two adjacent legs, the color-ordered amplitude A n is singular. (It is 
finite when the two collinear legs are not adjacent arguments.) This singular behavior has a universal form 
expressed by the tree-level factorization equation, 

AT e (l,...,a x %b x \...,n) k -^° Yl Spht^V*, 6^)^(1,..., (a + by,...,n), {22) 

ph. pol. cr 

dropping terms finite in the limit. In this equation, Split trcc is the usual tree splitting amplitude, and the 
notation 'a + b'' means k a + kb- The notation 'ph. pol.' indicates a sum over physical polarizations only. 
('Physical' here is in the sense of 'transverse', and their number may depend on the number of dimensions 
and on the variant of dimensional regularization employed.) This factorization is depicted schematically 
in fig. 1. It is characteristic of gauge theories that the splitting amplitude has a square-root singularity, 
Split <~ l/y/s a b, rather than a full inverse power of the two-particle invariant. 
At one loop, the analogous color decomposition to (2.1) is 

L"/2J+1 

A n ({k i ,X i ,a i })=g n J2nj E E G 

r„ ;c (cr) A^) c {a), (2.3) 

J c=l o-e5„/S„ ;c 

where \x\ is the largest integer less than or equal to x and nj is the number of particles of spin J. The 
leading color-structure factor, 

Gr„ ; i(l) = N c Tr(T ai ■■■T a ") , (2.4) 
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is just N c times the tree color factor, and the subleading color structures are given by 



Gr n;c (l) = Tr(T ai • • • T ac - X ) Tr (T a = • • -T a "). (2.5) 

S n is the set of all permutations of n objects, and S n[C is the subset leaving Gr„ ;c invariant. The de- 
composition (2.3) holds separately for different spins circulating around the loop. The usual normalization 
conventions take each massless spin-J particle to have two (helicity) states: gauge bosons, Weyl fermions, 
and complex scalars. (For internal particles in the fundamental (N c + N c ) representation, only the singlc- 
trace color structure (c = 1) would be present, and the corresponding color factor would be smaller by a 
factor of N c .) 

The subleading-color amplitudes A n - C> i are in fact not independent of the leading-color amplitude A n -i. 
Rather, they can be expressed as sums over permutations of the arguments of the latter [5]. (For amplitudes 
with external fermions, the basic objects are primitive amplitudes [6] rather than the leading-color one, but 
the same dependence of the subleading color amplitudes holds.) 




Figure 2. A schematic depiction of the collinear factorization of one-loop amplitudes. 

As a result, it suffices to examine the collinear limits of leading-color amplitudes. The collinear limits 
of the subleading-color then follow using this relation. The leading-color amplitudes obey the following 
factorization [5,11], 

^(1 n)-^ 



]T ( Split tre CT c (a A » , b x ") 4^(1, . . . , (a + by, . . . , n) 

ph. pol. a 

+ Split it° P («\ A^\{1, ...,(a + by,...,n) 
This factorization is depicted schematically in fig. 2. 
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(2.6) 



a 




a+b 



Figure 3. The defining equation for the one-loop splitting amplitude. 



This form was originally deduced from explicit calculations of higher-point amplitudes, but can also be 
proven more generally using the unitarity-based method. The latter proof also provides an explicit formula 
for the one- loop splitting amplitude [11], 

Split 1 _t° P (a Aa ^ Ab ) = 



E 



ph. pol. Ai,A2 



(2tt) 4 - 2£ I 2 



^ Splits + a + k)- A2 ,M)- Al ) 



(£ + k a + h) 2 



(2.7) 



x A\ rcc ((-£-a-b) x \a x ^b x \£ Xl ) 



The restriction to physical polarizations is important; it will give rise to transverse projection operators 
inside the loop. The defining equation is depicted graphically in fig. 3. 



3. A Sample Calculation 

As an example of a calculation using eqn. (2.7), we present the calculation of the splitting amplitudes 
Split 1_loop ( a + ; b~) in the conventional dimensional regularization (CDR) scheme, but restricted to physical 
external helicities. (The answer in the original 't Hooft-Veltman (HV) scheme would be identical. For a 
discussion of the differences of the various schemes, see refs. [15,16].) 

For this purpose, we need the three-gluon tree-level splitting amplitude, 

\/2 

Split trcc (P — > ab;z) = (-£„ • e b k b ■ e P +k b -e a e P -e b - k a -e b e a -e P ) . (3.1) 

Sab 

The variable z denotes the collinear momentum fraction, given in the collinear limit by 

q ' ka (3 2) 

q ■ (k a + h) 

where q is a reference null momentum which is not collinear with k a or k b . All momenta are taken to be 
outgoing, P + k a + k b = 0. The splitting amplitudes will in general have not only an explicit dependence on 
the argument z, but also an implicit one that arises from the dependence on the momenta and polarization 
vectors of legs a and b. Note that the order of the two arguments a and b is important, as the definition of 
z depends on it. We also need the four-point amplitude A 4 (l+, 2~, 3, 4), 
^ 4 (l + ,2-,3,4) = 
i 



2si 



(2-|^|l-)(2-|^|l-}-^- 4|l " )<2 " l ^ 3|2+) 



+ 



(12) (3.3) 

(2-|fe|l~)(l + |&fc|l~) 2g 3 - e 4 (2-|# 4 |l-) 2 
[12] «i 2 



where the notation indicates that the polarization vectors for the last two legs are kept generic, and in 
D = 4 — 2e dimensions. In this formula, s,j = (fcj + fcj) 2 , and following the notation of Xu, Zhang, and 
Chang [17], Ij*) and (j ± \ are massless two-component spinors carrying momentum kj; (j I) = and 
[jl] = are spinor products [18,17]. 

In the CDR scheme, there are 2 — 2e 'physical' polarizations, reflected in the identity 

Ect -a , li^v + k^q v 
=-9^ + ^ , (3.4) 

ph. pol. <T 

where k is the momentum of the gluon and q is the above-mentioned reference momentum. 

Using this identity, along with collinear identities such as q |a~) = — \fz{\. — z)2q ■ P, commuting 
gamma matrices appropriately, and dropping terms which are finite in the collinear limit, we find for the 
integrand of cqn. (2.7) a sum of terms, 

T (£) T (P -£) + Tl (i) - Tl (P-£)+ T 2 (£)] ^— 1__ , (3 . 5) 



where 



1 



T (p) = -=y/z{\-z)q-P{b-\t P \a-) 

, — ; p-P q- P (b~\dp\a~) (b~\^\a~) 

T 1 (p) = ^2z{l-z)q-P^ * ^ '^ P| lfl \ (3.6) 



TM . _ 2 ^(i ; ^^("-MI°-) 2 _ ^r\t- l«-> <»" 



The integrals we must perform then include various standard tensor two- and three-point integrals, 
along with integrals of the following form, 

f d 4 - 2e p 1 



Jl(Sab,z) 



(2tt) 4 -^ p\p-k a y{p-k a -k b y P -q> 



. I d^p 1 ,„ 



^ ^_ _. f d^p 



x) 2 (p- fc a - k b ) 2 p-q ' 

These integrals resemble those of axial gauge. We emphasize however that no explicit choice of gauge has 
been made above, and because all calculations leading to the cut expression are done on-shell (albeit in 4 — 2e 
dimensions), none is needed. Furthermore, the splitting amplitudes will be independent of q (as they must 
be). 

The integrals clearly scale like 1/r if we scale q — ► rq. They can depend only on q ■ k a and q ■ k b 
in addition to s ab . Other than the dot products of q with external momenta, however, s a b is the only 
dimensionful parameter. Furthermore, q ■ k a = z q ■ (k a + k b ) in the collinear limit, so that the first two 
integrals necessarily have the form 

Jl(Sab, Z) = 1 h( Z ) , 

(sab) q-(a + b) ^ 

•h{Sab, z) = — — f 2 (z) . 

(Sab) q- (a + b) 
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Consider first J\. To evaluate this integral, introduce Feynman parameters, and perform the loop 
integration, 



Ji{s ab ,z) = -12i J d 4 c6{l-^c<) J jt- 



1 



) 4 2e [-(1 - c 4 )p 2 + 2p ■ (c 2 k a + c 3 k a + c 3 k b - c 4 q) - c 3 s ab } 4 



-cic 3 s ab - 2c 4 q ■ (k a + k h ) (c 2 z + c 3 )] 



2+e 



(3.9) 

In order to solve for fi(z), we may set and q- (k a + k b )/ s ab to any convenient non-zero value. It will be 
convenient to choose it to be 1/2; at this point, 



1 



( 47r )2- £ 

Using the following change of variables, 



h{z) = -^(Sab) +£ Jl(Sab,z)\ 2q<ka+kb> 

r(2 + £) f d 4 c 5(1 - S 2 c 2 ) 



Sab 



(l-c 4 ) 2 



C1C3 + C2C4Z + C3C4] 



2+£ 



(3.10) 



ci=y(l-x), c 2 = (l-y)w, c 3 = (1 - y)(l - w) , c 4 = yx , 



(3.11) 



(with jacobian y(l — y)), we find 
T(2 + e) f 1 



Hz) 



r(2 + £ ) r 

(47r) 2 - £ J 



dxdydw [y(l - y)] 



-1-e 



(l-xyf 



[1 — w + iroz] 



2+e 



= ~2 Cr 



-r(l - e)r(l + e) z- x - e (l - z) £ - - + 2 Fi (e, e; 1 + e; z) 

z z 



1 , (l-z) e 



-r(i - e)r(i + e ) z- l - f -{\ - z y-- + 

z z 



+ — Li 2 (z) 

z 



(3.12) 



where 2-F1 is the Gauss hypergeometric function, Li 2 the dilogarithm, and the standard one-loop prefactor 
is, 



cr 



r(l + e)r^(l-e) 



(47r)2-er(l - 2e) ' 

The other integral in eqn. (3.7), J 2 , is manifestly z-independent, and one finds that 



(3.13) 



h{z) = — ^c r . 



(3.14) 



The last integral we need can be obtained by a Passarino-Veltman reduction [19,20], 



^3 { s abj z ) 



2 (s ab ) 1+e q-(a + b) 
With these results, we obtain immediately 
Split 1 - loop (P -> a+ &-;z) = 

- ^J~ b [zfiW + (1 - z)h(l z) 2/2] ( 



A(z) fc£ + -J- (2/ 2 - zA(z)) kl + -J- (A(z) - 2/2) e 

1 — z 1 — z 



(3.15) 



-S ab 



(1 ~ 2z) 



(6 I $ P \a ) + (fc„ - kb) ■ £ P 



(3.16) 



The corresponding tree-level result is, 



\/2s ah 



(1 ~ gf) 
^zJT^zj 



(b \ifp\a ) + {fa - fa) ■ £ P 



(3.17) 



4. General Form of Splitting Amplitudes 



We can also use the basic eqn. (2.7) for the one- loop splitting amplitudes in order to obtain a more 
general expression. We can repeat the steps in the previous section to obtain an expression for the gluon 
splitting amplitude with arbitrary external polarization vectors. In this case, we would use the general 
four-point amplitude rather than the expression (3.3), sum over polarizations crossing the cut using the 
identity (3.4), and perform the loop integrals. Doing so, we find (in the CDR scheme) the following result 
for the one-loop splitting amplitude, 



Split 



1-loop , 



P 



:b\: 



-Sab 



\zh{z) + (1 - z)h{\ -z)- 2/ 2 ] Split trco (P - ab;z) 



+ 



V2 S 2 b (l-2e)(3-2e) \-s ab 



h (fa - fa) ■ £P (SabSa ' £b - 2fa ■ E a k a • S b ) . 

(4.1) 

This calculation requires some tensor integrals beyond those presented in section 3; they are listed in ap- 
pendix I. As expected, the coefficient proportional to the tree-level structure has 1/e 2 and 1/e poles corre- 
sponding to soft and collinear singularities of the virtual gluons, as well as an ultraviolet singularity (the 
formula here is given before ultraviolet subtractions). The second polarization structure, which is new at 
loop level, is in contrast finite when e — > 0, as it must be. The two coefficients of the two polarization 
structures will make additional appearances later on, and so it will be convenient to define 



r 2 {z) 



2 \-s ab 

r 2 



[zf 1 (z) + (l-z)f 1 (l-z)-2f 2 ] 

,2 



(l-2e)(3-2e) \-s ab 



(4.2) 



The first term in eqn. (4.1) is in fact the same in all variants of dimensional regularization, but the second 
term is not. The corresponding results in the four-dimensional helicity scheme (FDH) [15] may be obtained 
by correcting for the difference (2 versus 2 — 2e) of helicity states with respect to the CDR scheme. This 
difference can also be expressed in terms of the contribution of an adjoint scalar inside the loop; for the 
general result, we find 



(1 - e5)e 2 



(l-2e)(l-e)(3-2e) 



-s ab 



h 



(4.3) 



where 



S = 



in the FDH scheme, 

1 in the CDR (or HV) scheme. 



We may write the separate contributions with particles of spin J circulating using the same polarization 
tensors, 



5Split[j] oop (P^a&;,j) 



r[ J] (z) Split troo (P a b; z) + r [ 2 J] (z) (fca ' ^ {s ab e a ■ e b ~ 2k b -e a k a - e b ) . 



ah 



(The formulae assume two helicity states, that is Weyl or Majorana fermions and complex scalars.) 

If we now take the external legs to be in four dimensions, we can write down expressions for the various 
helicity choices. For example, the following relations hold, 



Split 1 " 10015 (P+ -> a ± 6 T ; z) = Split trcc (P+ -> a ± b^; z) n(z), 

Split 1 " 10015 (P- -» a+ b+;z) = Split troo ( J p- -» a+ b+;z) ( n {z) - z(l - z) r 2 ), (4 5) 

,2 r 2- 



Split 1 - loop (P+ - a+ b+- z) = ^ z {l-z)\- ah] 



(ah)' 

The remaining helicity configurations can be obtained using parity. 

For completeness, we also record the contribution of a (complex) adjoint scalar circulating in the loop, 



rf ] = 0, 



J2- 



(l-2e)(l-e)(3-2e) \-s ab 

For the scalars in the fundamental representation of the gauge group, these functions must be multiplied by 
a factor of l/N. 

The basic cqn. (2.7) also applies to splitting amplitudes with fermions either circulating in the loop, 
or as external states. To compute the contributions to the gluon splitting amplitude proportional to the 
number of (massless) flavors n/, we need the tree- level g — > qq splitting amplitude (to the left of the cut), 
and the two-quark two-gluon amplitude (to the right of the cut), 

* Splits (P ^ab;z)= £ l / (gSl £ Splits - (< + a + % (-*),) {£ + ^ fc&)2 ^ 

x Af°(M-a-6) / ,a,6,£ / ) 

where the / subscript indicates fermionic legs. Evaluating this expression, we obtain for the coefficient of 
the two tensor structures, 

rl 1/2l =0, 

[1/2] = 6 2 (_£_Y . ( 4 - 8 ) 

7-2 " (1 - 2e)(l - e)(3 - 2e) v.- Sa J /2 ' 

The contributions proportional to the number of flavors is thus 

S Split 1 " 10015 (P ^ab;z) = ^ (fca ~ fcfc 2 ),£P ( Safc e a • £b - 2k b -e a k a - e b ) r^ /2] . (4.9) 

The l/N factor reflects the quarks being in the fundamental rather than the adjoint representation. 
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The one-loop g — > qq splitting amplitude requires, in addition to these tree-level splitting and four- 
point amplitudes, the four-quark amplitude. Here, there are two contributions which are independent of the 
number of flavors. One comes from gluons crossing the cut; the other, subleading in the number of colors, 
from quarks crossing the cut. The former is, 

5Split 1 - loop ( J P ^ a g bg;z) 



LC 



/j4-2ef a 



ph. pol. A1.A2 



(£ + k a + k b ) 2 



x Al ICC ((-£-a-b),a q ,bqJ). 



(4.10) 



The contribution from quarks crossing the cut is given by the following equation, 



6Split 1Aoop (P^a g bq;z)\ = 

I sc 



— y 

N 2 



ph. pol. A1.A2 



(27T) 



4-2e 



-Split tree (P^(* + a + 6) g (-*) ? ) 



(t + k a + k b f 



(4.11) 



x AT c ((-£- a- b) t ,a q >, bg,£ q ) 



where the subscript q indicates a quark, and q, an antiquark. Although the external quark and antiquark 
must necessarily have the same flavor, we make use of four-quark amplitudes with two different flavors in 
order to isolate the n/-independent contribution. We find for the two different contributions, 

<S Split 1 " 100 " (P-> a q bq;z) 



LC 



Split tICC (P^a q b q ;z)\r[ li (z)- 



1 



^S P lit 1 - l00P ( J P ^ flq bg-,Z) 



1 - 2e 

1; , tree/ 



(l-5e)e 2e(l-e)(l-£e) 
2(1 - 7) + (1 - 2e)(3- 2e) 



S a b 



(4.12) 



sc 



-Sab 



This splitting amplitude is again scheme-dependent; to compute the difference between the CDR and FDH 
schemes one can again either modify the gluon sums over polarizations or equivalently compute the contri- 
butions of (— e) complex internal scalars. In the latter method, the Yukawa coupling must be adjusted to 
obtain the 'same' coupling as for a gluon. The adjustment gives precisely the Yukawa coupling that would 
be obtained in an N = 4 supersymmetric gauge theory. The relevant scalar contribution, combining lcading- 
and sublcading-color parts, is 

*SpHtjoj oop (P->a,& 5; *) = 



Split^P-a,^;*) 



(l-2e)(3-2e) 2(1 - e)(l - 2e) iV 2 2(1 - e) 



Split troo (P^ a q bq;z) 



2e(l-e) 



1 



Sab 



h 



(4.13) 



.2(1 -e) (l-2e)(3-2e) iV 2 2(1 - e) / \ -s ab , 
In the first form, the last two terms arise from the Yukawa coupling. The second form matches the scheme- 
dependent parts shown in eqn. (4.12). In the above equations, the tree splitting amplitude is 

1 



Split" 00 (P -» a q bg-,z) 



V2s a , 



-UafpUb . 



(4.14) 
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Unlike the pure gluon case, for this splitting amplitude no other structure is possible when the fermions are 
kept on-shell. 

The contribution to the g — > qq splitting amplitude proportional to the number of flavors is, 



<f Split 1 " 10015 CLqbq-z) 



= 21 

v N 



x 



E / (^Sr j, ^ tICC (P + « + b) f (-/)/)- 



, .... (£ + k a + k b ) 2 

ph. pol. Al,A2 

xAT e {{-e-a-b)f,a q ,bg,£f) 



(4.15) 



We turn next to the fermion splitting amplitude, q — ► gg. All contributions here involve a tree-level 
g — > qg splitting amplitude to the left of the cut, and a two-quark two-gluon amplitude to the right; but 
there are both leading- and subleading-color ones. The leading-color term is, 

£ / (gSr j, Split- (P, - v + . + % (-0) - « 

ph. pol. Ai,A 2 

x 4 roo ((-£-a-%,a 9 ,6,£). 

(4.16) 

The result can be written in terms of the tree-level splitting amplitude, along with a new polarization tensor, 

Split troo (P (? -» a q b; z) = -i — u a j2* P u_ P , 

" V2 (4 ' 17) 
Split 1 " oop (P g — » a g 6; z) = Split troc (P 9 — » a q b; z)r\ c —u a ^bU-p k a ■ etr^ c . 

S ab 

Note that the momentum P is outgoing (in accord with our overall convention) , and hence U-p is the usual 
spinor for an inward-directed momentum. Given our conventions, helicity conservation means that the two 
external fermions at the ends of a given fermion line have opposite helicities: flipping the direction of a 
momentum also flips the sign of the helicity. 

The functions appearing in cqn. (4.17) are, 

^ = - *)/id - z) - ^-j:Kj ^ ( " 2 x 



2V (l-e)(l-2e)'7 \-s ab 



i , (4 - 18) 
„LC _ 1 (1 - 5e V ( A 4 ■ 



7-4 2(l-e)(l-2e) \-s ab ) h ' 

The subleading-color terms are given by a formula similar to (4.16), but with the fermion legs not 
adjacent in the four-point amplitude to the right of the cut, 

5Split 1 - loop (P g -» a q b;z)\ = 

j4-2c 



E / Ji SpHt trcc (^ - (/ + a + b) f (-*))- 



x< co ((-^- a -6) / ,a 9 ,£, 6). 
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(4.19) 



We may express the result of the sublcading-color contributions in a decomposition similar to that in 
eqn. (4.17), but multiplied by an overall factor of 1/N 2 . The coefficient functions are, 



2 V (l-e)(l-2e)'V \-s ab 



: , (4.20) 



1 (1 - 5ey ( fj, 
2(l-e)(l-2e) \s ab 

The difference between schemes again turns out to be expressible using contributions in which the gluon 
crossing the cut is replaced by a scalar; then for both the leading-color and subleading-color terms, these 
contributions are, 

3 " 2(l-e)(l-2e) {-s ab J /2 ' 

2 / 2 \ e l 4 - Zi 7 

JO] = £ f u 

4 2(l-e)(l-2e) 72 ' 

(The subleading-color contributions are of course again multiplied by a factor of 1/N 2 , not included here.) 
For specific helicity choices for the external legs, we then find, 

Split 1 " 10013 (P- -> a+ b~;z) = Split troo (F- - a+ b~; z)(r 3 + r 4 ) , 

(4.22) 

Split 1 " 1000 (P- -» a+ 6+; z) = Split trcc (p- -» a+ b+; z)(r 3 + zr 4 ). 

5. Consistency Checks 

The splitting amplitudes computed in the previous section may appear to be completely independent 
quantities. They are not. If we set n/ /N — > 1 and 1/N 2 — > —1, we obtain the corresponding quantities in a 
supersymmetric theory. These latter satisfy various identities descended from the supersymmetry Ward iden- 
tity [21]. Such identities thus provide a cross-check on the splitting amplitudes even in nonsupersymmetric 
theories such as QCD. 

Supersymmetry Ward identities have long been recognized [22,1] as providing useful relations between 
different amplitudes in gauge theories. For example, they relate amplitudes with n gluons to those with 
(n — 2) gluons and two gluinos. For certain choices of the external helicities, these identities, which hold 
to all orders in perturbation theory, allow us to solve for the two-gluino amplitudes in terms of the n-gluon 
ones. For example, 

Ar c (1^2-,3+,...,^,...,n+) = -MAr(l+,2-,3+,...r,...,n + ). (5.1) 

At tree level, the two-gluino color-ordered amplitudes are in fact identical to two-quark ones, so the identity 
links two different QCD color-ordered amplitudes. Beyond tree level, color-ordered QCD amplitudes are of 
course no longer identical to those in a supersymmetric theory. Supersymmetry identities arc nonetheless 
still useful in relating different contributions to QCD amplitudes [10]. 

The splitting amplitudes are in fact identical for all matterless supersymmetric gauge theories, indepen- 
dent of the number of supersymmetries. We may therefore pick an N = 1 supersymmetric gauge theory, 
with no matter content, to study the supersymmetry identities. 
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What relations do we expect for the splitting amplitudes? If we consider the n = 5 case of cqn. (5.1), 
with legs relabeled, we obtain the following relation between an all-gluon and a two-gluino thrcc-gluon 
amplitude, 

(3 1) A 5 (l+, 2+ 3~, 4+ 5-) = (3 5) A 5 (lp 2+ 3", 4+ 5 J) . (5.2) 

This relation holds to all orders in perturbation theory. Let us now examine the collinear limit k\ \\ k 5 . At 
tree-level, we then obtain the identity (P = — fci — fc 5 ) 

(31) J 4 trco (2 + ,3-,4 + ,(-P)-)Split troo (P+ -» 5- l+;z) = 

(5.3) 

(3 5) A trcc (2+, 3", 4+, (-P)-) Split tree (P+ -» 5g 1±; z) . 

(Note that the pure gluon amplitudes with a lone negative helicity vanish to all orders in the supersymmctric 
theory. The sum over the polarizations of the fused leg k\ + k$ thus reduces to a single term.) This identity 
in turn implies a relation between the g — > gg and g — > gg splitting amplitudes, 



Split troo (P+ -» a- q bp, z) = J- — - S P lit trcc (P+ -» a" &+;*), (5.4) 



where we have used k 5 — z(ki + k 5 ) and k\ = (1 — z)(fcx + /c 5 ) in the collinear limit to evaluate the ratio 
(31)/ (3 5). 

If we now repeat this exercise for the one-loop amplitudes, again taking the limit of eqn. (5.2), we obtain 
an identity of the same form for the loop splitting amplitudes, 



Split 1 - loop (P+ -> a~ bp, z) = 1— - Split 1 " 10013 (P+ ^a-b+;z). (5.5) 

The similarity is quite natural as the Ward identities hold to all orders in perturbation theory. 

One subtlety that arises in adapting the computations of the previous section to use in a supersymmetric 
theory is a difference in phase conventions. It is therefore more convenient to take the ratio of the tree-level 
and one-loop identities, 

Split 1 - 1000 (P+ -> a- g bp,z) = Split 1 -'°°P(P+ -> or b+; z) 

Split trcc (P+ ^ a- bp,z) ~ Split troo (P+ -> a-b+;z) ' ( ' ' 

as a test of our computation. This form of the identity is not dependent on the relative phase conventions 
for gluons and gluinos. 

In performing this check, we should of course use a supersymmetry-preserving regulator, so we pick the 
FDH scheme {6 = 0). 

To obtain the g — > gg splitting amplitude, start with the g — > qq amplitude (4.12,4.15), 



Split trco (P^ a q bg; z) 



1 { ) (l-2e)(3-2e) \-s ab ) h 



N 2 ) 

and set n f /N -> 1, l/N 2 -» -1, and 5 = 0: 



l-2e 2(1 -e) 



(5.7) 



Split 1 - 100 ^^ -^agbg-z) = Split trcc (P -» ag 6g; z) r[ 1] (z). (5.8) 
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In an N = 1 supersymmetric theory, the fermionic contribution to g — > gg cancels the r^ 1 ' term in the 
gluon contribution, and so we obtain, 

Split 1 " 10015 (P ^ab;z) = Split trcc (P -> ab:z)r l i ] (z) . (5.9) 

Both sides of eqn. (5.6) are thus equal to r^\z), and the identity is obeyed. 

The vanishing of + r^ 2 ^ is in fact also a result of a supersymmetry Ward identity. To see this, make 
use of the following Parke-Taylor equation [23] , 

A n (l+,2+,...,n-) =0, (5.10) 

which arc consequences of the Ward identities along with helicity conservation for fermionic amplitudes. 
Consider the collinear limit k\ \\ ky, of this equation at tree level. Using the above equation itself, the 
polarization sum reduces to a lone term (P = — (ki + £2)), 

4T((-P)-, 3+ . . . , n~) Split(P+ -> 1+ 2+; z) = . (5.11) 

Since the amplitude does not vanish, the splitting amplitude must. Repeating the above limit for loop 
amplitudes, we see that this result holds to all orders in perturbation theory. But in an N — 1 supersymmetric 
theory, this implies that 

Split(P+ - a+ b + ;z)\ {1] + Split(P+ a + b+;z)\ [1/2] = or r 2 1] = -r 2 1/2] . (5.12) 

If we consider theories with additional supersymmetries, or massless matter multiplcts in an TV = 1 theory, 
we can also relate the contribution of internal scalars to the fermionic contribution: 

Split(P+ -» a+b+;z)\ [1/2] + Split(P+ -» a + b+;z)\ m = or r 2 1/2] = -rf ] . (5.13) 

Inspection of cqns. (4.6) and (4.8) shows that the results presented herein indeed satisfy this relation. 

We can also derive a relation between the g — ► gg and the g — ► splitting amplitudes. To do so 
consider the collinear limit k 4 \\k 5 of eqn. (5.2) at tree-level (P = — k 4 — fc 5 ), 

(31) A trcc (l+,2+,3-,(-P)-)Split trcc (P+ ^4+5-;z) = 



(5.14) 

(3 5) A troo (l+ 2+, 3-, (-P)r) S P lit troo (P+ ^ 4+ 57; z) . 



Using the relation 



we obtain 



(3 1) 2+ 3- 4-) = (34) A(l+ 2+ 3",4r ) , (5.15) 



Split -» a+ ; z) = ^ Split (P+ -» a+ z). (5.16) 



1 

s/T^z' 

We have omitted the superscript 'tree' because repeating the limit for loop amplitudes we find that it holds 
for them as well. To avoid problems with phase conventions, it is again better to form ratios, 



Sp lit 1 - 1 °°P(P+ -» fl + b- g -z) = Split 1 - 1000 (P+ -> a+ b-;z) 

Split troo (P~+^a+^;z) ~ Split troo (P+^a+6-;z) ' ( ' ' 
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Setting l/N 2 



-1 in cqns. (4.17,4.18,4.20) we obtain, 

Split 1 " 10015 {P- a -» a g b; z) = Split troo (Fg -> a g b; z) r\ ] (z) 



(5.18) 



for the g — > 55 splitting amplitude, with r[ X '(z) defined in eqn. (4.2). As a result of setting l/N 2 — > —1, the 
function multiplying the structure 

u a ft b U- P k a -e b (5.19) 

(new at the one-loop level) vanishes. We are thus left with only one function needed to describe the one loop 
result for the g — > gg splitting amplitude, independent of the external helicities. We note that although one 
should use a supersymmetry-preserving regularisation scheme (i.e. FDH) it is remarkable that the difference 
between the results in the FDH and CDR schemes cancels out in the expression above. 

Comparing with eqn. (5.9), we see that this second supersymmetry identity is also satisfied. 

Given the observation that in supersymmetric theories, only one universal function multiplying the 
tree-level results appears in the one-loop splitting amplitudes, the two checks described above for a specific 
combination of helicities are sufficient to check the supersymmetry identities for all possible helicities. 

We may also compare our results with results previously obtained in the literature. In particular, we 
will compare with the result of Bern, Del Duca, and Schmidt [8] for the g — > gg splitting amplitude to all 
orders in e, and the results (to C(e )) for the g — > qq and q — > qg splitting amplitudes quoted in ref. [5]. For 
the comparison with ref. [8] it is convenient to use a slightly different representation for fi{z) than the one 
given in eqn. (3.10). Transforming the arguments of the hypergeometric function by the use of two identities, 
cqns. (9.131.1) and (9.132.2) of ref. [24] we may write 



-r(l-e)r(l + e) 



1 - z 



1-zl+e 



2 Fi(l,l + e;2 + e;— 



and 



(1 - z).h(l -z) = -\c T 2*1(1, -e; 1 - e; -?—). 

€ Z — 1 



(5.20) 



(5.21) 



We may now expand the hypergeometric functions in terms of poly logarithms, to obtain 



1-zl+e 
and thence, 



! 2^(1,1 + 6; 2 + 6;-^)- 2*1(1, -e; 1 - c; -Lj.) = 2 £ 

fe=l,3,5,. 



z-V 



z-V 



c k Li a 



-1, (5.22) 



,(^) < [- r < i -^ i -'(¥)' +2 t ^(^) 



fe=l,3,5, 



(5.23) 



where the polylogarithms are defined as follows [25], 

Lii(x) = — ln(l — x) . 



Li fc (a;) 



dt 



Lifc_i(t) (A; = 2,3,...) 



(5.24) 
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Comparing eqn. (5) of ref. [8] with eqn. (4.5), we see that comparing our results with theirs reduces to 
checking the following pair of equations, 

G n = rP(z), 



G f = -z{l-z) (r^+r-W) . 



(5.25) 



Comparing eqn. (5.23), and eqns. (4.3, 4.8) with the results presented in ref. [8], it is easy to see that the one- 
loop results shown in this paper indeed satisfy eqn. (5.25) and thus agree with the result^ given by Bern, Del 
Duca and Schmidt to all orders in the dimensional regularization parameter e. The contribution of an internal 
scalar is related to the difference between the FDH and CDR schemes in the gluon splitting amplitudes. Since 
our results agree with ref. [8] for both schemes, we implicitly agree on the scalar contribution as well. 

The representation of given in eqn. (5.23) is also useful for comparing the results for the g — > qq 
splitting amplitudes with the results quoted in the literature [5]. Using 



r(l- e )r(l + e) = l + - e 2 + 0(e 4 ), 



(5.26) 



we obtain from eqn. (5.23) 



r[ 1] (z) 



e 2 
cr 

c 2 



Sab 
2 



1 -Z, 



1 + eln(l - z) + eln(z) - — In ( 

2 z 



-e 2 + 0(e 3 ) 



-Sab 



1 - Z 



+ e 2 ln(z)ln(l-z)-— e 2 + 0(e 3 ) 



(5.27) 



for the expansion in e. Together with the expansion of the prefactors, 



2e(l 


-e)(l- 


-Se) 


(1- 


2e)(3- 


-26) 


1 


(1- 


Se)e 


1 - 2e 


2(1 
e 




(1- 


2e)(3- 


-26) 



= 1 



1 



10 2 

y ~ 3 

7 1 

2 + 2 



6)e 2 + 0(e 3 ), 



(5.28) 



= r+r 2 + 0(63), 



appearing in eqn. (4.12), (4.15) and (4.13), this leads immediately to the following expansions in e for the 
g —* qq splitting amplitudes, 



S Split 1 ~ loop (P — > a q bq\ z) 



= ^fc r Split troo (P -> a q bg\ z) 

Tlf iV 



2 1 



3 e \-s ab 



S Split 1 - 100 ^ - a q b q ; z) _ = -j^c r Split trcc (P - a q z) 



<5Split 1 - loop (P 



> a q b q ; z) 

Z(-Sab) 



SC 



LC 

6 

+ 



-Sab 



3 1 
2 7 



-Sab 



7 1 



0(e) 



c r Split trec (P^a 9 ^z 



(5.29) 



(1-Z)(-Sab) 
2 



A 4 



-Sa6 



13 1 

6 e 



A* 



1 



t Note that there is misprint in eqn. (9) of the journal version of ref. [8] due to a misplaced line break. One should replace 
rx (1 + e) by r(l + e). 
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For the contribution from virtual scalars we obtain, 



5 Split^ loop (P -» a q bg; z) = c r Split troo (P -» a 9 ^; z) 



11/ 



3 e \- s ah 



(5.30) 



if we consider only the first term in eqn. (4.13). (Recall that the two other terms in eqn. (4.13) are due to 
the coupling of the scalars to quarks.) These expansions in e agree with the results for the g — > qq splitting 
amplitudes presented in eqn. (B.12) of rcf. [5]. 

We close these comparisons with that for the q — > qg splitting amplitude. For this purpose, we must 



expand the functions r 3 
function zfi(z), 



LC.SC lc.SC • 



! '4 



in e. This can be easily done using the following expression for the 



e 2 Li 2 (z) + ln(z)ln(l-z)- 



+ 0(e), 



+ 0(e) 



e 2 Li 2 (l-z) 



(5.31) 



where we have used the dilogarithm identity 



U 2 (x) + Li 2 (l - x) = — - ln(x) ln(l - x) . 



(5.32) 



LC SC LC SC -i- 

The connection between the functions r 3 ' , r 4 ' used in this paper and the functions r s (q + ,a~), 
rs(q + ,a + ) introduced in eqn. (B.7) of rcf. [5] can be read off from eqn. (4.22): 



c r r s (q+, a-) = r^ c + r\ c + -jL ( r fC + r SC) 



c r r s (q\ o+) = r\ C + r 4 LC + ± (rf c + rf c ) - (1 - z) (r\ 
Using the expansion for zfi(z) we obtain 



„sc 



N 2 



(5.33) 



(1 - z)(-S a b) 



U 2 (z) 



cr 



N 2 



1 



e 2 \z(-Sab), 

-(! ~ *) + ^ 4 SC ) = cr(l + ^(1 -z) + O(e) . 
These results agree with eqns. (B.10) and (B.ll) of ref. [5]. 



£ \-Sab 



Li 2 (l-z) 



+ 0(e), ( 5 - 34 ) 



6. The One-Loop Antenna Function 



Gauge-theory amplitudes are also singular in the soft limit, when a gluon four-momentum vanishes. 
Color-ordered amplitudes display a factorization similar to that in eqn. (2.6), 



A^°°v(l,...,a,s x °,b,...,n) k " 



Soft trec ( a , S A % 6)4-i°° p (l,..., a ,6,...,n) 

+ Soft 1 - lo °P(a, S \ 6)^(1, . . . , a, 6, . . . , n) , 



(6.1) 
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in which the collinear splitting amplitude Split is replaced by the eikonal amplitude Soft. The cikonal limit 
also demonstrates one of the advantages of the looking at color-ordered amplitudes, namely that they allow 
a simple and straightforward factorization. The amplitude as a whole does not factorize simply, because 
color-charge factors get tangled up with the eikonal function. The eikonal function may be extracted to 
(9(e) from known results for five-gluon amplitudes. It has also been computed to all orders in e by Bern, 
Del Duca, and Schmidt [8]. 

In some applications, for example the construction of jet programs, one wishes to integrate over the 
singular regions. One can do this by integrating over the collinear and soft regions, and using in each the 
factorization appropriate to it. However, one must then introduce an artificial boundary between the two. 
It would be nicer to have a single, unifying function which summarizes the behavior in all singular regions. 

Catani and Seymour [26] wrote down just such a function, which appears in a dipole factorization 
formula. They presented a single function capturing the singular behavior of the squared matrix element 
in both the soft and collinear limits. One of us then gave a derivation of a single antenna factorization 
function [27] — so-called because it is within each color antenna that the two limits are combined — at 
the amplitude level. In this section, we will construct the one-loop analog of the antenna function and the 
antenna factorization formula. 

The antenna function presented in ref. [27] describes the factorization of an amplitude when a color- 
adjacent trio of legs a, 1, and b degenerates into two independent, hard momenta a and b, without a or b 
becoming soft. (This is the limit where A(a, 1, b)/s^ b — > 0, A being the symmetric Gram determinant of the 
(a, 1, b) system.) At tree level, the antenna function is given by the following formula, 

Ant(a,6 < — a, 1,6) = — e- h ■ e b V 3 (a, l,a) + — e a -e a V 3 (l,b,b) , (6.2) 
Sal s b i 

where V3 is the color-ordered gluon three-point vertex. In addition, we must specify a pair of reconstruction 
functions k- 1 — f~ t{k a , k\, k ) such that (k a ,k b ) — ► — (k a + ki,k ) sufficiently quickly when k\ \\ k a , and 
likewise (k a , k b ) — > — (k a , kb + fci) sufficiently quickly when ki || fc . These functions are not unique; different 
choices will lead to functional forms for the antenna function that differ by terms non-singular in any of 
the singular limits. We can also express the tree-level antenna function in terms of the tree- level splitting- 
amplitude as follows, 

Ant trcc (a, b <— a, 1, b) = -e b ■ e b Split troo (a - a 1; 1 - ^|) - e & ■ e a Split trcc (fe - 1 & ; , (6.3) 

where K 2 = s &b = (k a + ki + h) 2 . 

By analogy, we may now define the following function, 

L(a, b <— a, 1, b) = -e b ■ e b Split 1 " 10 ^ -> a 1; 1 - - e & ■ e a Split 1 " 100 ^ -> 1 6; ^) . (6.4) 

What are the limits of this function in the three singular regions, k a \\ k\, k b || fci, and k\ — > 0? 
In the first region, Split 1 ~ loop (6 — > lb; z) has at worst a logarithmic singularity as s a \ — > 0. Furthermore, 
kj, can be thought of as a reference momentum for the purposes of defining z, 

Sab $ab /n r \ 

z = ~ . (6.5) 

Sib + s ab S a l + s lb + s ab 
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The L function thus reduces to the ordinary one-loop splitting amplitude in this region. Similarly, it reduces 
to the splitting amplitude in the region fc \\ k\. 
In the soft region, k\ — > 0, 



Split trcc (a ->ol;z)-> k a ■e 1 s a -e a , 

Sal 

n 

Split trco (6^ lb-z) -y - — k b ■S 1 S b -St, 

Sib 



(6.6) 



while the second term in eqn. (4.1) is finite, and thus may be dropped. In taking the limit of L, we may 
set the relevant z to either zero or one as appropriate, except in expressions of the form z e . The leading 
behavior of L in the soft limit is thus given by 



'Sal 



i^a ' £l £a ' £d £o ' ££ 

•Sal 



A 4 



2 \ £ r 



-Sl6 



^2 

Sib 



fcfc ■ £l £ a ■ £a £6 ' £ 



6 ( ' 



using 



and 



|c r [-r(l - e)r(l + e) (1 - z) e - 1 + (1 - zf T{\ + e)r(l - e)] = 2/ 2 : 



z/i(z) z ~° -|c r r(i + c )r(i - e) z - f - 



(6.7) 

(6.8) 
(6.9) 



we obtain (K 2 = s a b in the limit) 



cr 



/ M 2 (-Sq b ) V W1 , > rn , /\/2 V2 \ 

7 77 r r(l + e)r(l - e) £ a • £a £{, • £H fc a ■ £1 h-Si) 

\ {Sal) {Sib) J Vs al Si6 / 



e 2 \(-Sai)(-sifc) 
c r / A* 2 (-s ah ) 



cr 



(-Sal)(-Slfc) 

A* 2 (-s afc ) 



T(l + e)r(l - e) £ a • £„ £ 6 • £ fe Soft trec (a, 1, 6) . 



■ £a • £a £fc • £ b Soft troo (a, 1, 6) . 



(6.10) 



e \(-Soi)(-si6)/ sm7re 
The latter form agrees with the result given in ref. [8]. 

The function we have defined is thus the appropriate generalization of the antenna function to one loop, 

Ant 1Aoop {a, 6 < — a, 1, b) = L(a, 6 < — a, 1, b) . (6.11) 

With it in hand, we may summarize the singular behavior of a leading-color one-loop color-ordered amplitude 
whenever the trio of momenta (a, 1, b) become degenerate without either a or b becoming soft as follows, 



Ant trcc (a, b <— a, 1, b) A^_°° P {1, . . . , -a(a, 1, 6), -6(6, 1, a), ... , n) 



ph. pol. (a, 8) 



(6.12) 



+ Ant 1 - loop (a,6^- a, 1,6) A^l, . . . , -a{a, 1, 6), -6(6, 1, a), . . . , n) ^ . 

As is the case at tree-level, this function contains terms which are not singular in any of the soft 
or collinear limits. Such terms lead to finite contributions when integrated over 'singular' regions of phase 
space. The inclusion or exclusion of such terms is arbitrary; in a physical observable they will cancel between 
'virtual+singular' contributions and 'real-emission' contributions. 
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7. Conclusions 



The cut-based method has been used extensively in recent years for computations of various gauge-theory 
amplitudes. We have shown here that it can also be used for simple, direct computations of the universal 
functions governing the collinear behavior of loop amplitudes, to all orders in the dimensional regulator e. Wc 
have also shown how to combine the different singular limits of real emission in one- loop amplitudes — soft 
and collinear into a single antenna function. A related function introduced by Catani and Seymour [26] for 
tree-level amplitudes, the dipole factorizing function, has already proven useful in a formalism for handling 
and cancelling infrared divergences in next-to- leading order calculations. The antenna function presented 
here should be similarly useful in next-to-next-to-leading order computations of jet processes in perturbative 



We thank Z. Bern, V. Del Duca, W. B. Kilgore, and C. R. Schmidt for an advance copy of a preprint [28] 
on another method for performing the calculations discussed here. 



In addition to the linear tensor integral, given in section 3, for the calculations in section 4, we also 
need the following integrals, 



QCD. 



Appendix I. Integrals 




Jg"(s ab ,z) = -i 




l J (2tt)4-2 £ pZ(p-k a y(p~k a -k b fp-q 



(1 - 2e) 2q ■ (a + b) 



1 s ab 




(1.1) 




-f 5ab (z) (kw + KK 



) + f 5bb (z)Kk v b 



Aq-{a + b) 



h bq (z) (fc 6 V + «) 



(ib 



2f5qq(z)q^q" + s a b.hg(z)g> 



+ 



4(g-(o + 6)) 
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where 



fbaa{z) 
faab(z) 
fbbb{z) 
haq{z) 
hbq{z) 

hqq{z) 
Ug{z) 



AW, 



2(1 - 2e) 

(l-e)z 



(l-2e)(l-z) 

2(l-2e)(l-z) 2 
e 

"(l-2e)(l-z) 



AW- ^"^I aW, 



(l-2e)(l-^- 

/l(z)+ (l-2 £ Kl-.) /2(z) - 
2(1 - ez) 



AW + 



2(l-2e)(l-z) 2 

z 



AW 



(l-2e)(l-z) 2 
(2-e-z) 



AW: 



AW , 



4(l-2e)(l-z) 



(l-2e)(l-^- 
/l(z) "2(l-2 C )(l-z) /2(z) - 



(1.2) 



For completeness, we also list the ordinary tensor bubble and triangle integrals of which we make use. 
These are of course independent of the collincar momentum fraction z, and may therefore be expressed using 
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1 2a 



1% 



1 



(2tt) 4 - 2£ p*(p-k a -k b ) 2 
e 1 



(1 - 2e) (-^fc) 1 



/2, 



P 



( 27r )4-2 £ p 2 (p _ ka _ kb) 2 

e 1 



2(1 - 2e) (-« o6 ) 

V (27T) 4 " 2 



■ /2 (fc„ + , 



2£ p2( p „ ^ _ /, b )2 

h \s ab g 



2(3-2e)(l-2e) (-* o6 )' 

j4-2e_ 



(2-e)(fc a + fc b )^(fc a + fc b ) 1 



Js = 



(-Sab) 



(27r) 4- 2e p 2 {p _ ka) 2 {p _ ka _ h) 2 



1 + c 



j4-2c 



P 



(2tt) 4 - 2 - p 2 { P ~k a y(p-k a -k b y 



(1-26) 

d 4 ~ 2e 



Sab) 
P 



-1-e 



(1-26) 



(-Sab)" 1 " 6 /2< , 



(27r) 4-2 £ p 2 {p _ ka) 2 {p _ ka _ kh) 2 

' ^ s afc / 2 . 9 ^ 



+ 



(-Sab) 1+£ 4(l-e)(l-2e) 
(2-6) 1 



2(1-26) (- 8ab ) 1+e h 1^ (fc " K + k " K) 



j4-2e 



P 



p H p v p P 



( 27r )4-2e p 2 {p _ ka) 2 {p _ ka _ kh) 2 

1 e(2-e) 



(-Sab) 1+£ 4(l-e)(l-2e)(3-2e) 

1 6 



sab/2 {g» v K + 9 w K + g vp K) 



+ 
+ 



(-s ah ) 1+e 4(l-2e)(3-2e) 
(2-e)(3-e) 1 



sab/2 WK + g^K + g^K) 



2(l-2e)(3-2e) (- Sab ) 

X f h^h u h p — 6 (h^h u h p A- h^h v h p 4- fFFl 
x ft a ft a _ l^a ^a ft b + ft a ft b ^a + K b ft a ^a 7 



2 - 6 



( h p h v h p 4- h p h v h p -4- h^ 1 h u h p \ 
\ K aH K b + K b ^"6 + ft b ""b^aJ 



—^—h p h v h p 

b 6 h 



(1.3) 
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